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Abstract 

Motion planning and control are key problems in a collection of robotic ap- 
plications including the design of autonomous agile vehicles and of minimalist 
manipulators. These problems can be accurately formalized within the lan- 
guage of affine connections and of geometric control theory. In this paper we 
overview recent results on kinematic controllability and on oscillatory controls. 
Furthermore, we discuss theoretical and practical open problems as well as we 
suggest control theoretical approaches to them. 

1 Motivating problems from a variety of robotic 
applications 

The research in Robotics is continuously exploring the design of novel, more reliable 
and agile systems that can provide more efficient tools in current applications such 
as factory automation systems, material handling, and autonomous robotic appli- 
cations, and can make possible their progressive use in areas such as medical and 
social assistance applications. 

Mobile Robotics, primarily motivated by the development of tasks in unreachable 
environments, is giving way to new generations of autonomous robots in its search for 
new and "better adapted" systems of locomotion. For example, traditional wheeled 
platforms have evolved into articulated devices endowed with various types of wheels 
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Figure 1: Underactuated robots appear in a variety of environments. From left to 
right, a planar vertical take-off and landing (PVTOL) aircraft model, a horizontal 
model of a blimp and the snakeboard. 

and suspension systems that maximize their traction and the robot's ability to move 
over rough terrain or even climb obstacles. The types of wheels that are being 
employed include passive and powered castors, ball- wheels or omni-directional wheels 
that allow a high accuracy in positioning and yet retain the versatility, flexibility and 
other properties of wheels. A rich and active literature includes (i) various vehicle 
designs |||, ||, |1 Egf, (ii) the automated guided vehicle "OmniMate" [§, (hi) 
the roller- walker ||l5|] and other dexterous systems [[l7| that change their internal 
shape and constraints in response to the required motion sequence, and (iv) the 
omni-directional platform in Jll| . 

Other types of remotely controlled autonomous vehicles that are increasingly be- 
ing employed in space, air and underwater applications include submersibles, blimps, 
helicopters, and other crafts. More often than not they rely on innovative ideas to 
affect their motion instead of on classic design ideas. For example, in underwater 
vehicle applications, innovative propulsion systems such as shape changes, internal 
masses, and momentum wheels are being investigated. Fault tolerance, agility, and 
maneuverability in low velocity regimes, as in the previous example systems, are 
some of the desired capabilities. 

A growing field in Mobile Robotics is that of biomimetics. The idea of this ap- 
proach is to obtain some of the robustness and adaptability that biological systems 
have refined through evolution. In particular, biomimetic locomotion studies the 
periodic movement patterns or gaits that biological systems undergo during loco- 
motion and then takes it as reference for the design of the mechanical counterpart. 
In other cases, the design of to robots without physical counterpart is inspired by 
similar principles. Robotic locomotion systems include the classic bipeds and multi- 
legged robots as well as swimming snake-like robots and flying robots. These systems 
find potential applications in harsh or hazardous environments, such as under deep 
or shallow water, on rough terrain (with stairs), along vertical walls or pipes and 
other environments difficult to access for wheeled robots. Specific examples in the 
literature include hyper-redundant robots Jll| 

snakes and roller racer models in |Bh. p7f , fish robots |23|, |25| , eel robots |37(] , and 
passive and hopping robots [jl8[ |36|, |43|. 

All this set of emerging robotic applications have special characteristics that pose 
new challenges in motion planning. Among them, we highlight: 
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Underactuation. This could be owned to a design choice: nowadays low weight 
and fewer actuators must perform the task of former more expensive systems. For 
example, consider a manufacturing environment where robotic devices perform ma- 
terial handling and manipulation tasks: automatic planning algorithms might be 
able to cope with failures without interrupting the manufacturing process. Another 
reason why these systems are underactuated is because of an unavoidable limited 
control authority: in some locomotion systems it is not possible to actuate all the 
directions of motion. For example, consider a robot operating in a hazardous or 
remote environment (e.g., aerospace or underwater), an important concern is its 
ability to operate faced with a component failure, since retrieval or repair is not 
always possible. 

Complex dynamics. In these control systems, the drift plays a key role. Dynamic 
effects must necessarily be taken into account, since kinematic models are no longer 
available in a wide range of current applications. Examples include lift and drag 
effects in underwater vehicles, the generation of momentum by means of the coupling 
of internal shape changes with the environment in the eel robot and the snakeboard, 
the dynamic stability properties of walking machines and nonholonomic wheeled 
platforms, etc. 

Current limitations of motion algorithms. Most of the work on motion plan- 
ning has relied on assumptions that are no longer valid in the present applications. 
For example, one of these is that (wheeled) robots are kinematic systems and, there- 
fore, controlled by velocity inputs. This type of models allows one to design a control 
to reach a desired point and then immediately stop by setting the inputs to zero. 
This is obviously not the case when dealing with complex dynamic models. 

Another common assumption is the one of fully actuation that allows to decouple 
the motion planning problem into path planning (computational geometry) and then 
tracking. For underactuated systems, this may be not possible because we may be 
obtaining motions in the path planning stage that the system can not perform in 
the tracking step because of its dynamic limitations. 

Furthermore, motion planning and optimization problems for these systems are 
nonlinear, non-convex problems with exponential complexity in the dimension of the 
model. These issues have become increasingly important due to the high dimension- 
ality of many current mechanical systems, including flexible structures, compliant 
manipulators and multibody systems undergoing reconfiguration in space. 

Benefits that would result from better motion planning algorithms for 
underactuated systems. From a practical perspective, there are at least two 
advantages to designing controllers for underactuated robotic manipulators and ve- 
hicles. First, a fully actuated system requires more control inputs than an under- 
actuated system, which means there will have to be more devices to generate the 
necessary forces. The additional controlling devices add to the cost and weight of 
the system. Finding a way to control an underactuated version of the system would 
improve the overall performance or reduce the cost. The second practical reason for 
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Figure 2: Vertical view of an omni-directional robotic platform with 6 degrees of 
freedom and 3 nonholonomic constraints JlJ, ^9). This device is capable of highly 
accurate positioning, high payloads, and high speed motion. In its fully actuated 
configuration, the robot is endowed with 6 motors at the three wheels and at the 
three joints (fit, fa, fa)- However, underactuated configurations can arise because of 
failures or intentional design. 

studying underactuated vehicles is that underactuation provides a backup control 
technique for a fully actuated system. If a fully actuated system is damaged and a 
controller for an underactuated system is available, then we may be able to recover 
gracefully from the failure. The underactuated controller may be able to salvage a 
system that would otherwise be uncontrollable. 

2 Mathematical unifying approach to the modeling 
of robotic systems 

Most of the robotic devices we have mentioned so far can be characterized by their 
special Lagrangian structure. They usually exhibit symmetries and their motion is 
constrained by the environment where they operate. In the following, we introduce 
a general modeling language for underactuated robotic systems. 

Let q = [q l ,-..,q n ) G Q be the configuration of the mechanical system and 
consider the control equations: 

if + T) k (q)qlq k = -M* ^ + fcj(<?)<f + l?(«)ui + • • • + Y^(q)u m , (1) 

where the summation convention is in place for the indices j, k that run from 1 to 
n, and 

(i) V : Q — ► M corresponds to potential energy, and kj(q)q J corresponds to damp- 
ing forces, 

(ii) {r*- fe : k = l,...,n} are n 3 Christoffel symbols, derived from M(q), the 
inertia matrix defining the kinetic energy, according to 

r fc = - M mk ( dMmj + ^™ - dM A 

11 2 \ dq* dqi dq m ) ' 
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where M mk is the (to, k) component of M 1 , and, 

(iii) {F a : a — 1, ... , to} are the m input co- vector fields, and {Y a = M~ 1 F a : a 
1 m} are the to input vector fields. 

Underactuated systems have fewer control actuators, m, than degrees of freedom 
n > to. Other limitations on the control signals u a might be present, e.g., actuators 
might have magnitude and rate limits, or they might only generate unilateral or 
binary signals (e.g., thrusters in satellites). 

The notion of affine connection provides a coordinate- free means of describing the 
dynamics of robotic systems. Given two vector fields X, Y, the covariant derivative 
of Y with respect to X is the third vector field Vx Y defined via 

BY" 1 

(VjeY)' = -qjX* +r i jh X*Y k . (2) 

The operator V is called the affine connection for the mechanical system in equa- 
tion (0). We write the Euler-Lagrange equations for a system subject to a conser- 
vative force Yq, a damping force k{q)q and m input forces as: 

in 

V 4 <? = Y a (q) + k(q)(q) + £ Y a {q)u a {t). (3) 

a=l 

Equation (||) is a coordinate-free version of equation (Q) . A crucial observation is the 
fact that systems subject nonholonomic constraints can also be modeled by means 
of affine connections. In the interest of brevity, we refer to jl0[ |3l[] for the exposition 
of this result and the explicit expression of the Christoffel symbols corresponding to 
the Lagrange-d'Alembert equations. 



The homogeneous structure of mechanical systems. 

The fundamental structure of the control system in equation (|^) is the polynomial 
dependence of the various vector fields on the velocity variable q. This structure 
affects the Lie bracket computations involving input and drift vector fields. The 
system (ph is written in first order differential equation form as 



d 


q 




~dt 







-T(q >q ) + Y (q)+k(q)(q) 



E 



U a {t) 



where T(q,q) is the vector with ith component T^ k (q)q> q k . Also, if x = (q,q) 



Z(x) 



-r(q,q) 





Y a (q) 



, and fc lift (x) = 



, Y^ix)^ 
the control system is rewritten as 

m 

x = Z{x) + Y™(x) + k Mt (x) + J2 Y™(x)u a (t) 





Hq)(q) 
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Let hi(q,q) be the set of scalar functions on R 2n which are arbitrary functions 
of q and homogeneous polynomials in {q 1 , . . . , q n } of degree i. Let Vi be the set 
of vector fields on IR 2n whose first n components belong to hi and whose second n 
components belong to hi+x- We note that these notions can also be defined on a 
general manifold, see J7|. 

We are now ready to present two simple ideas. First, all the previous vector 
fields are homogeneous polynomial vector fields for some specific value of i. Indeed, 
Z € V\, fc llft £ Pq, and F a llft € V-\. Second, since the Lie bracket between a vector 
field in Vi and a vector field in Vj belongs to Vi+j, any Lie bracket of the given 
relevant vector fields remains a homogeneous polynomial. In other words, the set of 
homogeneous vector fields is closed under the operation of Lie bracket. 

A consequence of this analysis is the definition of symmetric product of vector 
fields. We define the symmetric product between Yb and Y a as the vector field 
(Y a : Y b ) = (Y b : Y a ) given by 



(Y b : Y a y = (Y a : Y b y = -^Y b j + + T% (Y^ + Y*Yj 



Straightforward computations show that (Y a : Y b ) ult = [Y b m , [Z g ,Y^ lt ]}. This oper- 
ation plays a key role in nearly all the control problems associated with this class 
of systems: nonlinear controllability |l4|, ^|, optimal control p, |30|, dynamic feed- 



back linearization |Q, algorithms for motion planning and stabilization || Q 
etc 



A series expansion for the forced evolution starting from rest. 

The homogeneous structure of the mechanical control system (^), together with 
the symmetric product, set the basis to establish the following description of the 
evolution of the system trajectories starting with zero initial velocity || As- 
sume no potential or damping forces are present in the system. Let Y(q,t) — 
Y^a=x Ya{q)u a {t). Define recursively the vector fields Vk by 

Vi{q,t) = J o Y(q,s)ds, V k {q,t) = V k ^(q, s))ds. 

Then, the solution q(t) of equation satisfies 

q(t)=J2Vk(q(t),t), (4) 
fc=i 

where the series converges absolutely and uniformly in a neighborhood of go and over 
a fixed time interval t € [0, T]. This series expansion provides a means of describing 
the open-loop response of the system to any specific forcing. As we will see below, it 
plays a key role in several motion planning and control strategies for underactuated 
robots. 
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3 Existing results on planning for underactuated 
systems 

To design planning algorithms for underactuated robotic systems, we advocate an 
integrated approach based on modeling, system design, controllability analysis, dex- 
terity, manipulability, and singularities. These analysis concepts are fundamental 
for robust planning algorithms that do not solely rely on randomization or nonlin- 
ear programming. We do not suggest closed- form planning algorithms, rather we 
envision methods that combine the best features of formal analysis and of numerical 
algorithms. 

For reasons of space, we cannot present a detailed account of all existing results 
on motion planning for underactuated systems, and not even of the results obtained 
within the modeling approach proposed in Section |[ Therefore, we focus on two 
specific control methodologies for motion planning: decoupled planning algorithms 
for kinematically controllable systems, and approximate inversion algorithms based 
on oscillatory controls. 

Section |3.l| reviews decoupled planning algorithms that exploit certain differential 
geometric properties to reduce the complexity of the motion planning problem (still 
to be solved via numerical algorithms). The notion of kinematic controllability is 
extremely effective: trajectory planning decouples from being a problem on a 2n 
dimensional space to an n dimensional space. Furthermore, various state constraints 
can be neglected in the reduced space. For systems that are not kinematically 



controllable and that require oscillatory controls to locomote, Section 3.2 presents 
motion planning algorithms based on approximate inversion. Both design methods 
are closely related to recent results on nonlinear controllability [Bj3^], power series 
expansions two time-scales coordinate-free averaging p|, p5[, and nonlinear 

inversion algorithms || |34[ . 

The strengths of this methodology are as follows. Both methodologies provide 
solutions to the corresponding problems, i.e., point to point and trajectory planning. 
These analytic results do not rely on non-generic assumptions such as feedback lin- 
earization, nilpotency or flatness. The results are coordinate-free and hence widely 
applicable, e.g., to aerospace or underwater robotics settings. Both methods are con- 
sistent, complete and constructive (consistent planners recover the known solutions 
available for linear and nilpotent systems, and complete planners are guaranteed to 
find a local solution for any nonlincarly controllable system). 



3.1 Kinematic controllability for underactuated robots 

The following decoupling methodology was proposed in || to reduce the complex- 
ity of the motion planning problem. The method is constructive (only quadratic 
equations and no PDEs are involved) and physically intuitive. 

We consider as a motivating example a common pick-&-place manipulator: Fig. [| 
shows a vertical view of a three-revolute-joints device. We investigate planning 
schemes for this system when one of its three motors is either failed or missing. 
We present a decoupling idea to reduce the complexity of the problem: instead of 
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Figure 3: A three-revolute-joints device. It can be proven || that any two-actuator 
configuration of this system is kinematically controllable, i.e., one can always find 
two decoupling vector fields whose involutive closure is full-rank. 

searching for feasible trajectories of a dynamic system in R 6 , we show how it suffices 
to search for paths of a simpler, kinematic (i.e., driftless) system in R 3 . 

A curve 7 : [0, T] 1— > Q is a controlled solution to equation ([j]) if there exist inputs 
u a : [0, T] — ► R for which 7 solves ( |l|). To avoid the difficult task of characterizing all 
controlled solutions of the system (|l|), we focus on curves satisfying 7 = s(t)X("f), 
where A is a vector field on Q, and where the map s : [0, T] — ► [0, 1] is a "time- 
scaling" parameterization of 7. Such curves are called kinematic motions. 

We call V a decoupling vector field if all curves 7 satisfying 7 = s(t)V(j) for 
any time scaling s, are kinematic motions. This definition is useful for three rea- 
sons. First, V is decoupling if and only if V and VyU are linear combinations in 
{Yi, . . . , Y m }. Second, decoupling vector fields can be computed by solving (n — m) 
quadratic equations. Third, if enough decoupling vector fields, say Vi,...,V p , are 
available to satisfy the LARC, we call the system kinematically controllable. In the 
latter case, we can plan motions for the kinematic system q = 5Za=i w a(t)Va(q), and 
they will automatically be controlled curves for the original system (Q). 

3.2 Approximate inversion via small amplitude and oscilla- 
tory controls 

As in the previous section, the objective is to design motion planning and stabi- 
lization schemes for underactuated systems. We propose perturbation and inversion 
methods as widely applicable approaches to solve point to point and trajectory 
planning problems. Let us regard the flow map $ of equation Q over a finite time 
interval as a map from the input functions it; : [0, T] — > R to the target state x(T). 
The ideal algorithm for point-to-point planning computes an exact (right) inverse 
of $. Unfortunately, closed form expressions for $ _1 are available only as- 
suming non-generic differential geometric conditions (e.g., the system needs to be 
feedback linearizable, differentially flat, or nilpotent). Instead of aiming at "exact" 
solutions, we focus on computing an approximate inverse map using perturbation 
methods such as power series expansions and averaging theory. Although these tools 
are only approximate, the resulting algorithms are consistent and complete. 
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Oscillatory (high frequency, high amplitude) controls for trajectory plan- 
ning. 

We present the approach in three steps and refer to |35| for all the details. As first 
step, we present a recent coordinate-free averaging result. Let < e -C 1. Assume 
the control inputs are of the form 

1 ft 

U, — -Ui [ -, t 



e \e 



and assume they are T-periodic and zero-mean in the first variable. Define the 
averaged multinomial iterated integrals of Ui, . . . , u m as 



Uki,...,k m (t) = , , T , — ; / ( / ui(T,t)dr) ■ ■ ■[ I u m (r,i)dr) ds . 



fel / ps \ k„ 

{ / ii, (t *W 1 

fci! . . . k m \ J 

Let a, 6, c take value in {1, . . . , to}. Let k a (resp. k a b) denote the tuple (ki, . . . , k m ) 
with k c — S ca (resp. k c = 6 ca + S c b). Then, over a finite time q(t) = r(t) + 0(e), as 
e — » 0, where r(i) satisfies 

V,r = Y (r) + k(r)(r) + £ -^(t) - 27^ (t) (Y a : Y a )(r) (5) 



a<b 

As a second step, given z a (t), Zb c (t) arbitrary functions of time, we propose the 
following inversion procedure 

(i) take the functions ipN(afi){t) — V^N(a, b) cos(N(a, b) t), where (a, b) i— ► N(a,b) £ 
{1, . . . , N} is an enumeration of the pairs of integers (a, b), a < b. 

(ii) select the following controls in (||), 

u a (t, q) = v a (t, q) + ~W a ^-,t 
o-l 

c—1 c— a+1 

where v a (t, q) are still to be chosen. 

After computing the averaged iterated integrals of the oscillatory inputs w a (t/e,t), 
equation (|^) for the averaged system becomes 

m 

V,f - Y (r) + k(r)(r) + v a (t, r)Y a (r) 



a=l 



a=l a<b 
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Figure 4: Approximate trajectory tracking for an underactuated PVTOL model by 
means of oscillatory controls. The curve to be tracked is shown solid, and the various 
oscillating curves correspond to different values of the parameter e 



As a third and final step, assume that all the vector fields of the form (Y b : Yb) 
belong to span{y a }. Let a ab : Q — s- K be such that (Y a : Y a )(q) = J2 b <^ab{q)Y b {q), 
q £ Q. Select 



1 m / m \ 

~5>6a(?) 6-1+ £ (Z bc {t)f) 
6=1 V c=b+l / 



v a (t,q) = z a (t) 
Then, we have 

m mm 

5> a (t,r)y a (r) = ^zf(i)y„(r) +E^W^ : ^(r), 

a— 1 a— 1 a— 1 

which implies that eq. (||) takes the final form, 

rn 

Vrr = Y Q (r) + fe(r)(f) + ^ ^W^W + ^ z bc {t)(Y b : F c )(r) , 

a=l 6<c 

The averaged system now has more available control inputs than the original one. If 
the input distribution X = span{y a , (Y b : Y c )} is full rank, then the latter system is 
fully actuated (i.e., one control input is available for each degree of freedom). If the 
input distribution X contains a sufficient number of decoupling vector fields, then the 
system is kinematically controllable. In both cases, we have reduced the complexity 
of the motion planning problem. 

Remark 3.1 (Small amplitude algorithms based on series expansions) A re- 
lated approach to motion planning relies on small amplitude periodic forcing; see || 
The planning problem is solved by approximately inverting the series expan- 
sion describing the evolution of the control system (cf. Section H). This inversion 
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procedure is very similar to the one presented above. Based on it, one can estab- 
lish two simple primitives of motion to change and maintain velocity, while keeping 
track of the changes in the configuration. These primitives can then be used as the 
building blocks to design high-level motion algorithms that solve the point-to-point 
reconfiguration problem, the static interpolation problem and the local exponential 
stabilization problem. Fig. |5| shows two examples of the execution of these algo- 
rithms. 





Q 




.0" « 60 

v 2 (t) 



30 y 60 




Figure 5: Illustration of the motion planning algorithms via small amplitude periodic 
forcing for a simple planar body (left) and the blimp model (right). The errors in 
the final configuration are within the same order of magnitude of the input employed 



4 Open problems and possible approaches 

Immediate open questions arising from the above-presented results arc the following: 

Kinematic modeling and control. The current limitations are as follows: the 
design problem is now reduced to planning for a kinematic system with the addi- 
tional constraint of zero- velocity transitions between feasible motions. This addi- 
tional constraint leads to poor performance when coupled with current randomized 
planners |20| [24], ^| that switch frequently between the available motions. The 
zero-velocity switches also create problems for trajectory tracking controllers based 
on linearization, since the system loses linear controllability at zero- velocity. Finally, 
there is no notion of time-optimality for these kinematic motions and there is no 
way of dealing with systems where oscillatory inputs are needed for locomotion (see 
below for a discussion on this point). Motivated by this analysis, we identify the 
following open issues: 

(i) Develop a catalog of kinematically controllable systems, including planar ma- 
nipulators with revolute as well as prismatic joints, parallel manipulator, ma- 
nipulators in three dimensional space and in aerospace and underwater envi- 
ronments (accounting for the different dynamics in such settings). Some pre- 
liminary work in this direction can be found in Analyze and classify the 
singularities that these vector fields possess as a prerequisite step for planning 
purposes. 



11 



(ii) A (left) group action is a map tp : G x Q — ► Q such that i/)(e, q) = q, for all 
q £ Q, where e denotes the identity element in G, and ip(g, i/>{h, q)) — ip(gh, q), 
for all g,h € G, q S Q. Usually G C SE(n), and then the action describes 
a rigid displacement of some components of the robot. An interesting prob- 
lem would be to identify conditions under which decoupling vector fields can 
be found which are invariant under such group actions. When this is the 
case, motion plans can be designed exploiting established "inverse kinemat- 
ics" methods; see |}8[ Chapter 3]. This simplification eliminates the need for 
any numerical procedure if the robot moves in an un-obstructed environment, 
or further reduces the dimensionality and complexity of the resulting search 
problem in complex environments. 

(iii) To tackle the difficulties inherent with zero-velocity transitions, it would be 
appropriate to develop randomized planners which require as few switches be- 
tween decoupling vector fields as possible, and to develop trajectory tracking 
controllers for these systems able to adequately perform through the singular- 
ities. 

(iv) Another interesting idea would consist of switching between decoupling vector 
fields without stopping. In some sense, this is also related to the problem 
of developing transitions between relative equilibria. Relative equilibria are 
"steady trajectories" that the system admits as feasible solutions. This family 
of trajectories is of great interest in theory and applications as they provide a 
rich family of motions with the simplifying property of having constant body- 
fixed velocity. Relative equilibria for systems in three dimensional Euclidean 
space include straight lines, circles, and helices. Despite partial results, no 
method is currently available to design provably stable switching maneuvers 
from one relative equilibrium to another (or from one decoupling vector field to 
another without stopping). A necessary preliminary step toward this objective 
is to analyze the controllability properties of underactuated systems moving 
along a relative equilibrium or along a decoupling vector field. 

Small-amplitude and high-frequency controls. The current limitations are as 
follows. The implementation of the small amplitude approach requires the computa- 
tion and manipulation of high order tensors, and the approach has a limited region 
of convergence. The implementation of the oscillatory control approach presents 
difficulties in most physical settings because of the required high frequency, high 
amplitude inputs. Motivated by this analysis, we think that the following are inter- 
esting issues to explore: 

(i) For the small amplitude controls formulation, open questions include (a) inves- 
tigate tight estimates for the region of validity of the truncations (simulation 
studies suggest that there are better bounds than the conservative ones cur- 
rently available), (b) design base functions optimal with regards to region of 
convergence and appropriate cost criteria, (c) design inversion algorithms for 
systems that are not linearly controllable. The latter setting is equivalent to 
a non-definite quadratic programming problem, i.e., to the problem of finding 
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sufficient conditions for a vector-valued quadratic form to be surjective (see M 
for a discussion on this subject). 

(ii) For the oscillatory controls formulation, standing problems are (a) investi- 
gate the use of high-frequency bounded amplitude controls, (b) characterize 
approximate kinematic controllability and differential flatness via oscillations, 
(c) investigate physical settings in which oscillatory controls are natural control 
means, e.g., micro-electromechanical robots, (d) investigate extensions of this 
coordinate-free perturbation theory to discrete-time nonlinear systems, and to 
distributed parameter systems and partial differential equations. 

(iii) An ambitious program would consist of developing schemes that combine the 
proposed analytic methods with iterative numerical algorithms. One approach 
is via homotopy and level set methods (l} Q as schemes that overcome the 
limitations induced by the small parameter (small convergence region or high 
amplitude high frequency). A second direction is to use the planner based on 
small amplitude controls as a local planner inside a global search algorithm 
based on randomization; see 0] for some preliminary results on local/global 
planners. 
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